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ABSTRACT 
The most practical way to get a spatially broad and continuous measurements of 
the surface temperature in the data-sparse cryosphere is by satellite remote 
sensing. The uncertainties in satellite-derived LSTs must be understood to 
develop internally-consistent decade-scale land-surface temperature (LST) 
records needed for climate studies. In this work we assess satellite-derived 
“clear-sky” LST products from the Moderate Resolution Imaging 
Spectroradiometer (MODIS) and the Advanced Spaceborne Thermal Emission 
and Reflection Radiometer (ASTER), and LSTs derived from the Enhanced 
Thematic Mapper Plus (ETM+) over snow and ice on Greenland. When possible, 
we compare satellite-derived LSTs with in-situ air-temperature observations from 
Greenland Climate Network (GC-Net) automatic-weather stations (AWS). We 
find that MODIS, ASTER and ETM+ provide reliable and consistent LSTs under 
clear-sky conditions and relatively-flat terrain over snow and ice targets over a 
range of temperatures from -40 to 0°C. The satellite-derived LSTs agree within a 
relative RMS uncertainty of ~0.5°C. The good agreement among the LSTs 
derived from the various satellite instruments is especially notable since different 
spectral channels and different retrieval algorithms are used to calculate LST 
from the raw satellite data. The AWS record in-situ data at a “point” while the 
satellite instruments record data over an area varying in size from: 57 X 57 m 
(ETM+), 90 X 90 m (ASTER), or to 1 X 1 km (MODIS). Surface topography and 


other factors contribute to variability of LST within a pixel, thus the AWS 
measurements may not be representative of the LST of the pixel. Without more 
information on the local spatial patterns of LST, the AWS LST cannot be 
considered valid ground truth for the satellite measurements, with RMS 
uncertainty ~2°C. Despite the relatively large AWS-derived uncertainty, we find 
LST data are characterized by high accuracy but have uncertain absolute 
precision. 


INTRODUCTION AND BACKGROUND 


Recently, there has been a great deal of attention on the increasing temperatures 
in the Arctic due to observed and projected warming (IPCC, 2007) and the 
positive ice-albedo feedback. Climate models predict continued Arctic warming 
but they differ in their predictions of the extent, rate and magnitude of the 
temperature increases. The most practical way to get a spatially broad and 
continuous measurement of surface temperature of the cryosphere is through 
satellite remote sensing. The uncertainties in satellite-derived land-surface 
temperatures (LSTs) must be understood to develop internally-consistent 
decade-scale LST records needed for climate studies. The level of confidence in 
the accuracy of LST records may be established through comparison with 
independent observations. LST can also be referred to as ice-surface 
temperature (IST) but in this work, we adopt the conventional LST terminology. 


Analysis of the Advanced Very High Resolution Radiometer (AVHRR)-derived 
surface temperature of the Antarctic and Greenland ice sheets, ice caps and 
smaller glaciers provides a method to evaluate melt and climate trends over 
about the last three decades. These measurements can be used as a surrogate 
for, and in addition to, air-temperature records (Box, 2002) that are obtained from 
relatively few and scattered in-situ observations (Stroeve and Steffen, 1998; 
Wang and Key, 2003 and 2005a and b; Comiso, 2006). Satellite-derived LST 
data are also useful for validation of climate models and as input to data- 
assimilation models. North of 60°, clear-sky surface temperature has increased 
over about the last two decades. Based on analysis of AVHRR data acquired 
above the Arctic Circle, Comiso (2006) shows an increase of 0.72+0.10°C per 
decade from 1981 — 2005; Wang and Key (2005b) show an increase of 
0.57+0.02°C per decade using the AVHRR Polar Pathfinder dataset from 1982 — 
1999. Warming has not been uniform, and there are some areas of cooling 
during that same period, particularly over large areas of Siberia. 


Thermal infrared (TIR) sensors facilitate LST and melt-condition monitoring over 
extensive areas (Key and Haefliger, 1992; Rees, 1993; Key et al., 1997; Stroeve 
and Steffen, 1998; Comiso, 2006; and Hall et al., 2004, 2006 & 2008), especially 
when used with complementary satellite-derived passive- and active-microwave 
data (e.g., see Nghiem et al., 2001). Yet a rigorous assessment of the random 
and systematic errors in the various satellite TIR measurements of surface 
temperature has not been established for snow and ice targets. LST errors 
derived from satellites come from a variety of sources, including: effects of very 
thin, unmasked clouds and atmospheric water vapor, and imprecise accounting 
for emissivity variations. In addition, calibration inconsistencies in the longer- 
term satellite record are important, especially when the same basic instrument is 
used on different satellites over a period of many years. 


NASA's EOS provides space-borne instruments suitable for measuring snow/ice 
surface temperature using TIR sensors. The present study focuses on TIR 


surface-temperature measurements from standard products derived from two 
satellite-borne instruments, namely: the Moderate-Resolution Imaging 
Spectroradiometer (MODIS) and the Advanced Thermal Emission and Reflection 
Radiometer (ASTER), both on the EOS Terra satellite. We also include LSTs 
derived from the Enhanced Thematic Mapper Plus (ETM+) on the Landsat-7 
satellite. These sensors have different spatial resolutions and spectral channels; 
in addition, different algorithms are used to derive surface temperature. Though 
MODIS and ASTER LSTs have been compared in many previous studies [see 
for example, Jacob et al. (2004), Schmugge and Ogawa (2006) and Mihalcea et 
al. (2007)], the studies have not focused on the cryosphere. To our knowledge, 
no study has been published in the peer-reviewed literature that compares 
multiple snow and ice LSTs derived from the various Earth Observing System 
(EOS) standard products. 


One way to assess the accuracy of remotely-sensed LSTs is to compare the 
values with in-situ LSTs. Various factors make that a difficult task including the 
fact that the in-situ observations are point measurements while the satellite- 
derived observations represent LSTs from a much larger area. Also, in-situ 
observations from weather stations, and in particular automatic weather stations 
(AWS) are generally acquired at some height above the surface so that the 
measured air temperature must be extrapolated to a surface temperature, or 
LST. That calculation is also affected by meteorological factors such as wind 
speed, but for higher wind-speed conditions, near-surface air temperature is very 
close to the surface temperature. 


Standard LST products, available from the Land Processes Distributed Active 
Archive at the EROS Data Center in Sioux Falls, South Dakota, are available 
from both the MODIS and ASTER instruments. These products were developed 
and validated using a specific algorithm that is community-endorsed for use with 
a specific sensor. Though other products may also be available, the standard 
products are usually maintained by the investigator who developed the algorithm. 
A great deal of documentation and quality-control information is available for the 
products from the product developers and the EROS Data Center. The 
validation presented in the present paper is useful for assessing the errors of 
these commonly-used products. Though there is no standard LST product from 
the ETM+ instrument, we use the formulation of Barsi et al. (2005) to calculate 
the LST. 


In this paper, we assess LSTs from MODIS and ASTER standard products 
(MOD11_L2 and MOD11A1, and AST08), and LSTs from ETM+ for snow and ice 
targets on the Greenland Ice Sheet, and ice caps located outside the margins of 
the ice sheet. We use in-situ air-temperature data from Greenland Climate 
Network (GC-Net) AWS (Steffen et al., 1996; Steffen and Box, 2001) (Figure 1) 
to compare with satellite-derived LSTs. It is not the intent of this work to derive 
more accurate LSTs nor to modify or “improve” existing algorithms, rather, we 
are comparing satellite-derived LSTs relative to each other and, relative to in-situ 


measurements, to improve our understanding of the uncertainties in, and the 
limitations of, these data products. These commonly-used products are 
increasingly considered for uses in climate-data records (CDRs), thus the 
uncertainties must be established. 


THEORETICAL BASIS FOR REMOTE SENSING IN THE MID- AND 
THERMAL-INFRARED TO RECOVER SNOW AND ICE SURFACE 
TEMPERATURES 


Mid-IR and TIR-remote sensing is based on measurement of the radiance 
emitted from the surface and modified by the atmosphere, with a goal to recover 
the kinetic temperature and emissivity from the measured at-sensor radiance. 
The at-sensor radiance (Ls) for a given wavelength (A) in the mid or TIR, 
excluding any reflected solar contribution in the mid IR, can be written as: 


Lsa= [ €a Lova(T) +(1- €2) Leta] Ta + Lema [1] 


where: 

€, = surface emissivity at wavelength 2. 

Lma(T) = spectral radiance from a blackbody at surface temperature T. 

Lissa = Sky radiance (spectral downwelling radiance incident upon the surface 
from the atmosphere.) 

tT, = transmittance (spectral atmospheric transmission.) 

Lama = path radiance (spectral upwelling radiance from atmospheric emission 
and scattering that reaches the sensor.) 


The effect of z, is to reduce the amount of ground-emitted radiance measured at 
the SENSOF; Lam ,adds a component unrelated to the ground and L,,,, serves to 


reduce the spectral contrast. Any surface radiance, temperature, or emissivity 
retrieval must compensate for atmospheric effects. 


The surface temperature (7) is not an intrinsic property of the surface; it varies 
with external factors such as meteorological conditions. Emissivity /s an intrinsic 
property of the surface and is independent of the temperature. See Hook et al. 
(2007) for further discussion. 


Surface emissivity 


The surface emissivity is defined as the ratio of the actual radiance emitted by a 
given surface to that emitted by a black body at the same kinetic temperature. 
Salisbury et al. (1994) show that snow emissivity departs significantly from black 
body behavior in the 8 — 14 um part of the spectrum. Emissivity of snow varies 
with liquid water content and snow grain size especially at larger grain sizes 


(Salisbury et al., 1994; Wald, 1994; Hori et al., 2006). Grain sizes on the surface 
of the Greenland Ice Sheet do not vary much due to redistribution of surface 
snow by the wind, allowing a constant grain size (0.3 — 0.5 mm) to be an 
appropriate assumption for non-melting snow according to Stroeve et al. (1996). 
[Other, special circumstances, can cause the surface snow-grain size to increase 
when, for example, surface hoar forms with centimeter-size snow grains 
(Shuman and Alley, 1993)]. 


Emissivity is also sensitive to sensor viewing angle (Dozier and Warren, 1982; 
Salisbury et al., 1994; Hori et al., 2006). Fine- and medium-grained snow 
exhibits only a slight angular dependence, but the angular dependence increases 
for coarse-grained snow, especially at the larger viewing angles (e.g., 75°). Bare 
ice exhibits the largest angular dependence, showing a decrease in emissivity 
with increasing viewing angle (Hori et al., 2006). Also see Dozier and Warren 
(1982) for modeled values of hemispherically-averaged emissivity curves of snow 
for grain radii from 50 — 1000 um, and for viewing angles of 0 to 75°, and Hori et 
al. (2006) who reviewed and discussed those results, and Dash et al. (2002) for 
further discussion of land-surface temperature and emissivity estimation. 


Emissivity of snow and ice features is far less variable than the range 
encountered in a wide variety of land surfaces, leading to greater potential 
accuracies in retrieved LST over snow and ice even when the emissivity is not 
known precisely. To obtain snow-surface temperatures to an accuracy of 0.1 K, 
the emissivity must be known to within 0.1% (Stroeve et al., 1996). 


Spectral libraries are available that provide emissivities based on laboratory 
measurements (Shafer, 1971; Salisbury et al., 1994; the ASTER spectral library 
http://speclib.jpl.nasa.gov/); and the University of California at Santa Barbara 
emissivity library (http:/Awww.icess.ucsb.edu/modis/EMIS/html/em.html)). Shafer 
(1971) measured the emissivity of 185 snow samples at temperatures from -1 to 
-18°C in the 8-14 um wavelength range, obtaining values ranging from 0.966 to 
0.989, with fresh snow having the lowest value (0.975), and crusted snow having 
the highest mean-measured emissivity (0.985), under clear skies. The average 
emissivity of all of his samples was 0.978, though more-recent works tend to 
report higher emissivity values (see, for example, Wan et al., 2002; Hori et al., 
2006). 


SATELLITE INSTRUMENTS 

MODIS 

The Terra satellite was launched in December 1999, and the Aqua satellite was 
launched in May 2002, both having MODIS instruments in their payload. The 


Terra and Aqua MODIS instruments are effectively the same, although band 31 
(11 um) on the Aqua MODIS saturates at temperatures about 60K below those 


for Terra (~400K). (In addition, the Aqua MODIS near-infrared band 6 is non- 
functional.) 


The MODIS instrument acquires data in 36 spectral channels (Barnes et 
al.,1998; Salomonson et al.,1989). MODIS scans + 55° from nadir and provides 
daytime and nighttime imaging of any point on the Earth every 1-2 days. 
Channels 31 and 32 centered on 11.03 and 12.02 um, respectively, are used to 
produce the LST products that are used in this study (Table 1). MODIS IR 
radiances are quantized in 12 bits and have a ground resolution of 1 km at nadir. 
MODIS IR radiances are calibrated with a cold space view and full aperture black 
body viewed before and after each Earth view. 


ASTER 


The ASTER instrument on the Terra satellite acquires radiances in 14 spectral 
channels as part of three sensor sub-systems. Channels 10-14 of the TIR 
subsystem (between 8.15 and 11.65 um) are used to produce the product that is 
used in this study (Table 1). The TIR data are quantized in 12 bits, have a 
spatial resolution of 90 m at nadir and are calibrated with a full-aperture 
temperature-controlled black body. The black body is viewed before and after 
each data acquisition and periodically heated to obtain calibration data at 
different temperatures. The TIR subsystem utilizes a whiskoroom scanner that 
can be pointed + 8.55° across track allowing any coordinate on Earth to be 
imaged as frequently as every 16 days. 


ETM+ 


The Landsat-7 ETM+ is also part of the EOS suite of satellites. Landsat-7 was 
launched in April 1999, and has eight bands ranging from 0.45 to 12.5 um; the 
ground resolution ranges from 15 m in the panchromatic band to 57 m in the 
single TIR band (Table 1), from which LST may be calculated. The ETM+ scan- 
line corrector failed in June 2003 rendering subsequent ETM+ data useless for 
the present work. 


Ice and Cloud Elevation Satellite (ICESat) 


The ICESat satellite, launched in January 2003, carries a single instrument, the 
Geoscience Laser Altimeter System (GLAS), which was designed to use laser 
pulses at 1064 and 532 nm to map the surface topography of Earth precisely and 
to sound the atmosphere. The primary I|CESat mission is to provide detailed 
surface-elevation data of the ice sheets, to improve ice-sheet mass-balance 
assessments and enable volumetric change monitoring (Zwally et al., 2002). For 
this work, we use ICESat-derived elevations along transects across two ice caps 
(Flade Isblink and North Ice Cap) outside the margins of the Greenland Ice Sheet 
to study variation in LST as a function of elevation. 


INSTRUMENT SPECIFIC LAND-SURFACE TEMPERATURE (LST) 
ALGORITHMS, PRODUCTS AND GROUND MEASUREMENTS 


The methodology used to derive LST varies from MODIS, ASTER and ETM+, as 
the different methods take advantage of the unique characteristics of each 
sensor. The accuracy of both the split-window and single-channel methods, 
discussed below, relies on knowledge of the surface emissivity and atmospheric 
temperature and humidity. Because snow/ice emissivity does not vary widely, 
these conventional methods produce useful snow/ice LSTs at the upper limit of 
accuracy. 


MODIS products 


We use the LST Collection 4 (or Version 4 (V4)) product suite (MOD11 and 
MYD11) developed by Wan et al. (2002); also see Wan (2008). (MOD refers to a 
Terra MODIS product and MYD refers to an Aqua MODIS product.) We use only 
Terra MODIS products because there are more opportunities for comparison with 
other data since the Terra MODIS has been operating longer than the Aqua 
MODIS. Furthermore, Terra has an overpass time (14:30-16:30 UTC) for 
Greenland within 2 hours of local solar noon. Also, MODIS and ASTER are both 
flown on the Terra spacecraft, so it is often easy to find MODIS scenes when 
ASTER scenes are available. 


The Terra MODIS 1-km resolution LST algorithm produces a swath product 
(MOD11_L2) that is based on the generalized split-window technique (Wan and 
Dozier, 1996) using MODIS bands 31 and 32 (Table 1). To determine the actual 
LST from an instrument that has two IR channels, one must correct for 
absorption and reemission of radiation by atmospheric gases, predominately 
water vapor. The “split-window” method is widely used to achieve some 
correction for atmospheric effects because the measured temperature difference 
between the two IR channels, in a “split-window” algorithm, is proportional to the 
amount of water vapor in the atmosphere. To compute the LST using 
MOD11_L2, the emissivity must be prescribed. For bands 31 and 32, the 
emissivities used in the algorithm to compute LST over the Greenland Ice Sheet, 
are 0.993 (for band 31) and 0.990 (for band 32), and do not vary seasonally. The 
daily 1-km LST product (MOD11A1) is developed from the swath product. 
MOD11A‘1 is provided on a sinusoidal grid. 


Wan and Dozier (1996) calculated coefficients in the generalized split-window 
LST algorithm by interpolation on a set of multidimensional look-up tables 
obtained by linear regression of MODIS simulation data from radiative-transfer 
calculations over a wide range of surface and atmospheric conditions. The 
coefficients depend on column water vapor and surface air-temperature values 
from the MODIS atmosphere-profile standard product, MOD07 (Seemann et al., 


2006). The MODIS Atmospheric Profile product consists of several parameters: 
total-ozone burden, atmospheric stability, temperature and moisture profiles, and 
atmospheric water vapor. All of these parameters are produced day and night for 
Level 2 at a resolution of 5x5 1-km pixels when at least nine FOVs are cloud free. 
The following MODIS products are input to the MOD11_L2 algorithm: sensor 
radiance (MOD021km), geolocation (MOD03), cloud mask (MOD35_L2), 
atmospheric temperature and water vapor (MODO7_L2), land-cover (MOD12Q1) 
and snow-cover (MOD10_L2). 


The swath-based LST product (MOD11_L2) was used in the present work when 
it was necessary to know the exact time that the data were acquired to 
synchronize image data and to coordinate with ground and other satellite 
observations. The daily LST product (MOD11A1) was used in the transects 
across the two ice caps where knowledge of exact time-of-day was not needed 
because the ice cap surfaces were at the melting point. 


ASTER products 


ASTER produces one emissivity product (AST05) and a separate surface 
temperature product (AST08). The emissivity product contains five bands 
corresponding to the surface emissivity at the wavelengths of the five TIR bands. 
The data are in swath (scene) format. A more detailed description of the ASTER 
instrument, products and their application is given in Yamaguchi et al. (1998). 
For this work, we used the AST08 products downloaded between 25 July 2006 
and 3 October 2007. 


The ASTER algorithm to derive LST and emissivity uses a radiative transfer 
model to correct the at-sensor radiance to surface radiance followed by an 
emissivity model to separate the surface radiance into temperature and 
emissivity. The ASTER LST algorithm requires atmospheric profiles from either 
satellite sounding or conventional weather balloon soundings made twice-daily, 
globally, for the atmospheric model (Price, 1983; Susskind et al. 1984; Chedin et 
al. 1985) and an emissivity model that is typically based on laboratory and field 
measurements (Kealy and Hook, 1993; Matsunaga, 1994). In this approach, a 
temperature-independent index is calculated showing the range, that is 
minimum-maximum difference (MMD) in emissivity with wavelength, but not the 
absolute value. A calibration curve is then derived using laboratory data that 
allows the observed variation (or MMD) to be used to obtain the unknown 
absolute emissivity value. 


Once the unknown emissivity value is obtained, Equation 1 can be used to obtain 
the temperature and emissivities for the other bands. Rocks and soils have a 
large MMD whereas vegetation and snow have a small MMD. This is because 
within the TIR range, rocks and soils show strong emissivity variation, whereas 
the emissivity of snow and ice is more consistent, varying between about 0.97 
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and 1.0. The calibration-curve approach is used to recover LST and emissivity 
from ASTER (Gillespie et al., 1998 and1999). 


ETM+ derived LST 


A separate atmospheric correction must be done, and the target emissivity must 
be prescribed, as in Barsi et al. (2003 and 2005), to derive LST from the single 
TIR band on the ETM+. Ancillary data are needed to adjust the top-of- 
atmosphere (TOA) brightness temperature to a surface skin temperature, or LST. 
An atmospheric-correction tool was developed by Barsi et al. (2005) for the 
Landsat-5 Thematic Mapper (TM) and Landsat-7 ETM+ thermal channels, and is 
accessible through http://atmcorr.gsfc.nasa.gov. Date, time and location of the 
Landsat scene are input by the user. The atmospheric-correction parameter 
calculator uses National Centers for Environmental Prediction (NCEP) modeled 
global atmospheric profiles as input data, along with MODTRAN and other 
integration algorithms, to derive site-specific atmospheric transmission and 
upward and downward radiances for the given Landsat scene. Using the 
atmospheric-correction parameter output, a series of calculations is performed to 
derive a surface kinetic temperature for each pixel in the Landsat scene. The 
first calculation involves converting the digital number, DN, to the TOA spectral 
radiance, R (see Markham and Barker, 1985). Then the formulation of Barsi et 
al. (2005) is used to perform the atmospheric correction, deriving the surface- 
leaving radiance, R. Finally, conversion of radiance to LST can be made using a 
sensor-specific approximation to Planck’s function. 


T = K2/ log [(K1/R) + 1] [2] 


Where T is surface temperature in Kelvins, and K2 and K1 are calibration 
constants in Kelvins, and in W m* sr‘ um", respectively (Markham and Barker, 
1985). 


Greenland Climate Network (GC-Net) Automatic Weather Stations (AWS) 


The GC-Net AWS network (Figure 1) has collected weather and climate 
information on the Greenland Ice Sheet since 1995 (Steffen et al., 1996; Steffen 
and Box, 2001) [http://cires.colorado.edu/science/groups/steffen/qcnet/]. Each 
AWS is equipped with instruments to sample micro-meteorological variables 
used to obtain surface energy and mass budget closure. Air temperature is 
measured at two heights within 5 m of the surface with 1.2 m separation using 
thermocouples and thermistors in passively-ventilated solar radiation shields. 
Samples are obtained at 15 or 60 second intervals. Hourly average data are 
transmitted via a satellite link throughout the year. Other measured parameters 
include: wind speed and direction, air humidity, air pressure, surface 
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accumulation rate, and surface radiation balance in broadband shortwave and 
TIR wavelengths. 


Continuous near-surface air temperature and LST records have been acquired 
for more than two years at the Summit AWS site using a Vaisala CS-500 Type-E 
TIR thermocouple sensor to record air temperature and a Vaisala Type-T 
thermocouple TIR sensor that was placed 10 cm above the surface to record the 
surface temperature as one of the GC-Net AWS stations 
(http://cires.colorado.edu/science/groups/steffen/gcnet/) assuming a snow 
emissivity of 0.99. An absolute accuracy for LST was found to be +0.2°C (r? = 
0.98) at the Summit site (Figure 2) where a mean warm bias of 1.0°C was 
evident in 16,123 hourly samples. The surface temperatures were derived from 
two positions on an AWS mast with actively-ventilated air temperature sensors. 
Field and laboratory calibrations of the thermocouple devices are generally 
conducted annually and the resulting corrections are typically <3% (Steffen et al., 
1996; Steffen and Box, 2001). 


METHODOLOGY AND RESULTS 


To assess the accuracy of MODIS and ASTER LST standard-data products, and 
LSTs derived from the ETM+, first we compared satellite-derived LSTs along 
transects across ice caps on Greenland that were experiencing surface melt 
across the entire ice cap. We then compared MODIS, ASTER and ETM+ LSTs 
of the ice sheet with AWS-derived LSTs when possible at a variety of surface 
temperatures. Finally, we compared MODIS and ASTER, and MODIS and ETM+ 
LSTs in the absence of in-situ observations, at surface temperatures ranging 
from about -16 to 0°C. 


Transects 


Transects over two ice caps in Greenland (located outside the margins of the ice 
sheet), Flade Isblink and North Ice Cap, were studied using MODIS and ASTER 
LST products, and ICESat altimetry data (Figures 3a and 3b). On the image 
dates (3 July 2001 and 23 June 2004), LSTs do not vary with elevation (derived 
from ICESat), thus the entire ice surface is approximately at the melting point 
(Figure 3a). The mean LST of Flade Isblink (not shown) on 3 July 2001 is 
0.84+0.70°C (n = 1,564,335), as determined from ETM+ (where the value 0.98 
was input as the emissivity), and 0.47+0.26°C (n = 21,843), as determined from 
MODIS (MOD11A1), an insignificant difference of only 0.37+0.44°C. MODIS and 
ASTER LSTs agree over North Ice Cap (23 June 2004) within ~+0.8°C (Figure 
3b). Agreement is also about ~+0.8°C between the satellite-derived LSTs and 
the surface if we assume that the surface is at 0°C. 


Comparisons of satellite-derived LSTs with GC-Net AWS-derived 
measurements 
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First, scenes were selected manually from the various instruments that were 
acquired on the same day, and at the same time (when possible). We searched 
for relatively-flat areas as close to nadir as possible, to minimize angular effects. 
Pixels from each scene were selected from homogeneous parts of the scene as 
determined from visual (textural) analysis, and in areas with no obvious clouds or 
fog. The LST value from one MODIS pixel was compared to mean LSTs of the 
same area from the ASTER and ETM+ scenes. 


We used the LST value from the MODIS pixel (MOD11_L2) corresponding to the 
location of the AWS and the time of the AWS measurement. Initially, we 
examined LST image pairs over the Greenland Ice Sheet and in-situ data from 
nine stations in the AWS network (Figure 4) for a total of 11 cases. The number 
of ETM+ pixels used to calculate the mean LST increases with distance from 
nadir, corresponding with the increasing size of the MODIS pixel away from 
nadir. Ten of the cases involved comparing MODIS and ETM+ LSTs with an 
AWS-derived LST, and one case (2 July 2001) involved comparing ASTER and 
ETM+ LSTs with the corresponding AWS-derived LST. (The MODIS LST 
product was not available on 2 July 2001.) In that one case where near- 
simultaneous ASTER and ETM+ were compared with the AWS-derived LST from 
Summit station, we used just one ASTER and one ETM+ pixel, the pixels that 
corresponded most closely with the location of the AWS. Since it was not always 
possible to get the times of the satellite-scene acquisitions and AWS acquisition 
of LST within an hour of each other, the AWS-derived air temperatures were 
interpolated between acquisition times to provide a more-accurate comparison 
with the time of acquisition of the MODIS LSTs (see Table 2). 


For the cases shown in Table 2, the LSTs derived from different satellite 
instruments are in better agreement than they are with the AWS LSTs. Also, the 
satellite-derived LSTs agree in sign in most cases (Figure 4). 


MOD11_L2 LST comparison with AWS observations. A collection of 121,506 
hourly-average cases from 15 GC-Net AWS (Figure 1) over the 1998-2001 
period are used to evaluate the near-surface (below ~5 m) temperature gradient 
(AT/Az) to model LST. AT/Az is measured using T,-—T,, where T, is farther 
from the surface than 7. Az is the instrument spacing: z,—z,, is equal to 1.2 m 
in most cases. The magnitude of AT/Az , that is, the near-surface thermal 
stratification varies with surface radiation balance and micro-meteorological 
conditions. Under clear-sky conditions, the strongest modulators of AT/Az are 
downward solar irradiance (SL) and wind speed. When SJ is small, surface 
emitted IR radiation exceeds absorbed solar irradiance and the surface cools 
more than the air above the surface leading to a temperature inversion, that is, 
AT/Az > 0. Wind speed facilitates turbulent mixing of heat downward under 
inversion cases where it is lost by upward IR emission. During air temperature 
lapse conditions, the surface has been heated, usually by SJ such that turbulent 
mixing moves heat away from the surface leading to negative AT/Az. Strong 
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winds can neutralize AT/Az even when surface solar heating is large. A clear 
signal (r=-0.987) of decreasing temperature inversion strength with increasing 
downward solar irradiance is evident. Below 20 W m®, on average, AT/Az= 
+0.33 K m'. On average, inversion becomes lapse for downward solar 
irradiance exceeding 350 W m*. The gradient in AT/Az with SJ is -0.90 K m' 
kW''. Thus, the lapse rate at 800 W m* = -0.40 Km". 


Figure 5 illustrates the comparison between MODIS-derived LST and AWS air 
temperature with and without cases associated with large AWS air temperature 
uncertainty, that is, when downward solar irradiance exceeded 240 W m®, wind 
speed was < 4ms', and air temperatures above the melting point. When 
neutral thermal stratification (that is, AT/Az| < 0.1 K m'') conditions prevail, 
nominal height air temperature may be used as an LST surrogate because air 
temperatures and LSTs are equivalent within the specified accuracy (~0.1 K) of 
the AWS air temperature sensors. 


Estimating LST by extrapolation during strong AT/Az cases is not useful for the 
MOD11_L2 accuracy assessment. It is only under small absolute A7/Az that an 
insignificant MODIS LST error is evident. 


It was not possible to compare significant numbers of ASTER and ETM+ -derived 
LSTs with corresponding AWS-derived LSTs (as was done with MODIS) due to 
our inability to locate an adequate number of cloud-free ASTER and ETM+ 
scenes coincident with automatic weather stations. 


Comparison of Satellite-Derived LSTs 


ASTER and ETM+ LSTs were compared with MODIS LSTs in the absence of in- 
situ measurements. These results provide additional information on how well the 
independently-derived LSTs agree over a range of ice-surface temperatures 
(also see Figures 3a and b for additional cases when LSTs were near 0°C), even 
though different instruments, algorithms and channels are used to derive the 
LSTs. Points were selected from a range of elevations on the ice sheet and in 
different seasons, covering a range of surface temperatures from about -16 to 
0°C. 


ASTER- and MODIS-derived LSTs agree, with r* = 0.97 and RMSE = 0.54°C (n = 
65). The mean bias = 0.22°C, with the ASTER LSTs being slightly higher than 
the MODIS LSTs (Figure 6). 


Results of the MODIS and ETM+ comparisons are shown in Figures 7a and b. 
Two different emissivity values were input to the ETM+ algorithm and were 
selected based on the literature survey. Using an emissivity value of 0.98 as an 
input to the ETM+ algorithm, r* = 1.0, (RMSE = 0.34°C, mean bias = 0.07°C (n = 
80), with ETM+ values being slightly higher than MODIS values. When an 
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emissivity = 0.99 was input to the ETM+ algorithm, r? = 0.99, RMSE = 0.56°C 
and mean bias = -0.43°C (n = 80), with ETM+ values being slightly lower than 
MODIS values. Though the correlation between ETM+ and MODIS LSTs is 
approximately the same and very high (r=1.0), the mean bias and RMSE are 
both lower when using an emissivity value of 0.98 (as compared to when 
emissivity = 0.99 was used). Implicit in that, however, is the assumption that the 
MODIS LSTs are the “truth,” but we have no way of knowing that, given the data 
that we have available. 


ADDITIONAL CONSIDERATIONS 
Clouds 


The satellite-derived LSTs represent “clear-sky” surface temperatures. ‘Clear 
sky’ is determined from the MODIS cloud mask (Ackerman et al., 2002) for the 
MODIS products used in this work. Though we studied each scene visually to 
confirm that the selected points were not cloud contaminated, undetected thin 
clouds and a varying atmospheric water vapor column may also influence the 
accuracy of satellite LST retrievals. Incorporating information from other 
sensors, such as the Atmospheric Infrared Sounder (AIRS) instrument could help 
constrain error from unidentified thin cloud and atmospheric water vapor radiative 
effects. A cloud-clearing algorithm that uses texture and time-variance might 
also be useful such as that reported by Shi et al. (2007) using the MODIS and 
Multi-angle Imaging Radiometer (MISR), both flown on the Terra satellite. 


Over Greenland, the effect of water vapor absorption on the brightness 
temperature in MODIS band 31 (Table 1) is estimated to be within +1K in the 
winter, and greater (2-3K) in the summer and in southern Greenland, especially 
around the coast. When the atmospheric temperature is higher than the LST, the 
atmospheric effect may increase the brightness temperature (Z. Wan, written 
communication, 2008). 


Angular effects 


Angular effects that may influence the accuracy of the MODIS LST retrievals 
were not studied in this work, nor were effects associated with rugged terrain. 
Based on modeling by Dozier and Warren (1982), Key and Haefliger (1992) 
show the emissivity variations by scan angle for the AVHRR channel 4 (11pm). 
For example, the emissivity varies from 0.9988 at a scan angle of 0°, to 0.9955 at 
a scan angle of 50°. However, the effect of scan angle is greater with coarse- 
grained snow and bare ice (Hori et al., 2006). In addition, ice-sheet topography 
can influence the emissivity and thus the accuracy of the LST retrieval. Though 
we attempted to select only relatively-flat areas, areas on the ice sheet that 
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appear flat from above can actually have significant local topography, that is, with 
a topographic undulation amplitude of >20 m. 


How representative are the AWS measurements? 


The AWS “point” measurements may fail to represent the surface temperature 
over the area of the 1 km x 1km (MODIS), 90 m x 90 m (ASTER) or 57 mx 57m 
(ETM+) pixels, especially in instances where wind speed is low and surface slope 
heterogeneity is relatively large. One author (K. Steffen) has made LST 
measurements over a 1-km area at Swiss Camp on the Greenland Ice Sheet and 
found surface-temperature variability on the order of ~1°C. The temperature 
variability was attributed to surface roughness (typically 10-15 cm in amplitude 
and several meters in wavelength) and some shadowing. Additional surface LST 
observations need be made in a local array to evaluate spatial variability in LST, 
over areas corresponding with satellite pixels. Without more information on the 
local spatial patterns of LST, the AWS LST cannot be considered absolute 
ground truth for the satellite measurements. 


SUMMARY AND CONCLUSIONS 


MODIS and ASTER products, and ETM+ derived LST provide reliable and 
consistent LSTs under clear-sky conditions and relatively-flat terrain over snow 
and ice targets over a broad temperature range of -40 to 0°C. The satellite- 
derived LSTs agree better with each other than they do with the GC-Net AWS- 
derived LSTs even though different algorithms are used with each instrument to 
calculate LST. 


The agreement between MODIS (MOD11_L2) and ASTER (AST08) LSTs is very 
good over the Greenland Ice Sheet, with r* = 0.97 and RMSE = 0.54°C (n = 65). 
The agreement is even better between MODIS and ETM+ especially when an 
emissivity = 0.98 was input to the algorithm to calculate LST, with r* = 1.0 and 
RMSE = 0.34°C (n = 80). 


We have demonstrated the difficulty in validating LST measurements from 
satellites using in-situ measurements even over a relatively homogeneous land 
surface such as an ice sheet. Careful selection of AWS data was necessary to 
obtain a reliable set of in-situ data to compare with the satellite data. The 
average bias of satellite-derived and in-situ LSTs (-0.3 °C) is well within the RMS 
error of 2.1°C (n = 48). Without more information on the local spatial patterns of 
LST, the AWS air-temperature-derived LSTs cannot be considered absolute 
ground truth for the satellite measurements because the AWS record data at one 
“point,” while the satellite instruments record data in a pixel varying in size from 
57 X 57 m (ETM+), 90 X 90 m (ASTER) or 1 X 1 km (MODIS). Surface 
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topography and other factors contribute to variability of LST within a pixel, thus 
the AWS measurements are often not representative of the entire pixel. 


It is not possible from this work to determine which LST product or measurement 
provides the highest accuracy, or which LST algorithm is the most robust. An 
extensive network of surface-temperature measurements, coincident with 
satellite overpasses, is required to answer those questions. 


In the future, the National Polar-orbiting Operational Environmental Satellite 
System (NPOESS) Visible/Infrared Imager/Radiometer Suite (VIIRS) will be 
launched with IR channels and an LST product will be produced. Other 
instruments may also be launched that will provide LST measurements or 
products. To facilitate the production of data records that will be developed using 
multiple instruments and perhaps even climate-data records (CDRs) over snow 
and ice targets, uncertainties in LST retrieval must be established. The use of 
internally-consistent products will permit determination of “clear-sky” surface- 
temperature trends (if any) to be detected and monitored in the polar regions for 
climate and other studies. As such, when the data are applied in climatological 
analyses, trends should exceed our calculated RMS errors to be considered 
‘significant’ within the observational capability. 
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TABLES 


1. Spectral range and spatial resolution of TIR bands on MODIS, ASTER and 
ETM+ that are used to derive the LSTs discussed in this paper. Earth Science 
Data Type (ESDT) refers to the name of the MODIS and ASTER standard 
products. 


2. Data corresponding with Figures 4, AWS and ASTER, MODIS and ETM+ LST 
comparison. There are 10 AWS-MODIS comparisons and 11 AWS-ETM+ 
comparisons; one ASTER and one ETM+ pixel each, are used for the one 
ASTER and ETM+ comparison. AWS acquisition times were interpolated to 
match MODIS acquisition times. 
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FIGURES 


1. Greenland Climate Network (GC-Net) automatic weather station (AWS) 
location map of Greenland (from Steffen and Box, 2001). 


2. Hourly air temperatures vs. surface temperatures in 2000 and 2001 for 
Summit Station, Greenland. 


3a. Map view of two specific ICESat altimetry profiles across the Flade Isblink 
ice cap in northeastern Greenland. The approximate center coordinates of Flade 
Isblink ice cap are 81° 34’ N, 15° 0’ W. Elevations from ICESat Track 0045 were 
acquired on 02/25/04 and elevations from ICESat Track 0248 were acquired on 
03/10/04 during Laser 2B operations period. The surface elevations derived from 
these ICESat tracks are shown in the above right graphs, and the black vertical 
bars show the approximate start and end of the ice surface. Plotted with the 
elevations are surface temperatures derived from the MOD11A1 product for 3 
July 2001. The entire ice cap surface is close to the melting point, so the surface 
temperature does not vary with elevation. 


3b. ASTER true color (left) and land surface temperature (LST) (right) images of 
North Ice Cap on 23 June 2004 are shown. ASTER and MODIS surface 
temperatures derived from along the transect shown on the ASTER LST image 
are plotted. The transect starts in the north and ends in the south. The 
approximate center coordinates of North Ice Cap are 76° 55’ N, 68° 00’ W. 


4. Comparisons of satellite-derived and in-situ AWS LSTs. MODIS-derived and 
in-situ AWS LSTs are compared for 10 cases; ETM+ and AWS-derived LSTs are 
compared for 11 cases; and one case shows a comparison of ASTER and in-situ 
AWS LSTs. Four MODIS and four corresponding ETM+ scenes, and one 
ASTER scene were used to develop the graph. 


5. Illustration of the comparison between MODIS-derived LST and AWS-derived 
air temperature in cases where the AWS air temperature should equal the LST. 
Excluded are cases when downward solar irradiance exceeded 240 W m* while 
wind speed was <4 ms". Also excluded are cases when AWS air temperatures 
were above the melting point. The solid line is the 1:1 line. The dashed line is 
the best-fit line. 


6. MODIS — ASTER comparison. Seven MODIS and seven corresponding 
ASTER scenes were used to develop the plot using three scenes from the 
Greenland Ice Sheet. 65 cases are plotted. The solid line is the 1:1 line. The 
dashed line is the best-fit line. 


7a. MODIS and ETM+ comparison, where emissivity = 0.98 in the ETM+ LST 
calculation. The solid line is the 1:1 line. The dashed line is the best-fit line. 


7b. MODIS and ETM+ comparison, where emissivity = 0.99 in the ETM+ LST 
calculation. The solid line is the 1:1 line. The dashed line is the best-fit line. 
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Table 1. Spectral range and spatial resolution of TIR bands on MODIS, ASTER 
and ETM+ that are used to derive the LSTs discussed in this paper. Earth 
Science Data Type (ESDT) refers to the name of the MODIS and ASTER 
standard products. 


Band number and |Pixel size 
(spectral range) in 


ee 


31 (10.780 - 11.28) |1 X 1 km 
32 (11.770 - 12.27) 

10 (8.125 —8.475) [90 X 90m 
11 (8.475 — 8.825 
12 (8.925 — 9.275 
13 (10.25 — 10.95 
14 (10.95 — 11.65 
6 (10.31 — 12.36) 7X57m 


No wm was Na 


7/7/2001 
8/1/2001 
6/29/2000 


8/3/2001 
6/8/2000 


7/2/2001 
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Table 2. Data corresponding with Figure 4 - AWS and ASTER, MODIS and ETM+ LST 
comparison. There are 10 AWS-MODIS comparisons and 11 AWS-ETM+ comparisons; 
only one ASTER and one ETM+ pixel each, are used for the ASTER and ETM+ 


comparison. AWS acquisition times were interpolated to match MODIS acquisition 


times. 
Distance 

AWS | Tor No. of ait: rasa ETM+ MODIS from 
AWS Station | LST (©) Pixels 2 LST (©) Avg a Acquisition Acquisition Nadir 

(°C) Avg. Avg. Averaged LST (°C) Time (UTC) Time (UTC) ese 
Crawford Pt 1 1306 14:50 15:25 79.05 
Crawford Pt 4 1258 14:43 15:15 30.28 
DYE-2 1271 14:34 15:15 93.54 
NASA-E -10.78 | -10.16 1251 -10.05 14:13 14:50 
NASA-SE -15.00 | -11.98 1219 -12.75 14:07 14:40 
NASA-U 1282 15:28 16:05 54.72 
NGRIP -10.08 1257 -10.65 15:03 15:40 | 0.00__| 
Saddle -14.42 | -14.37 1290 -15.37 14:19 14:50 64.58 
South Dome 1265 14:10 14:50 21.15 
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Figure 3b 
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Figure 5 
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Figure 7b 
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